Loewner driving functions for off-critical percolation clusters 
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We numerically study the Loewner driving function Ut of a site percolation cluster boundary 
on the triangular lattice for p < pc- It is found that Ut shows a drifted random walk with a 
finite crossover time. Within this crossover time, the averaged driving function (Ut) shows a scaling 
behavior —{pc—p) 4'"+^'/^" with a superdiffusive fluctuation whereas, beyond the crossover time, the 
driving function Ut undergoes a normal diffusion with Hurst exponent 1/2 but with the drift velocity 
proportional to (pc ~ pY , where v = A/Z is the critical exponent for two-dimensional percolation 
correlation length. The crossover time diverges as [pc ~ p)~^'^ as p ^ p^. 

PACS numbers: 05.40.-a,64.60.ah 



Introduction. Loewner evolution has recently drawn 
much attention in physics community because of the 
development of Schramm-Loewner evolution (SLE) 
which has provided us a new tool for the study of two- 
dimensional (2-d) continuous phase transition. The basic 
device of SLE is a conformal mapping that transforms a 
motion along a non-intersecting 2-d curve into another 
motion along the real axis; Loewner driving function is 
a real function that represents this transformed motion. 
It turns out that, for a certain class of stochastic and 
conformally invariant curves in 2-d, the driving func- 
tion shows Brownian motion in one dimension. What 
makes SLE especially remarkable is that it gives us a 
method that describes all the geometrical properties of 
the curves through a single parameter of the Brown- 
ian motion, namely, the diffusion constant. The class 
of curves includes the self-avoiding walk, the uniform- 
spanning trees, the loop-erased random walk, and bound- 
aries of critical clusters in various 2-d lattice models in 
physics such as percolation, Ising model, 0{n) loop mod- 
els, and Potts models (see [1| for review). 

Being inspired by mathematically oriented develop- 
ment, people start using SLE formalism to test the con- 
formal invariance by calculating Loewner driving func- 
tions obtained from 2-d curves in a number of physical 
systems, such as vorticity clusters and temperature iso- 
lines in turbulence Q , domain walls in 2-d spin glass [3| , 
isoheight lines on growing solid surface 5] , nodal domains 
of chaotic maps ^], etc. 

In study of physical systems, it is important to ask how 
driving functions may look like when the system departs 
from the critical point, because in real life there are a 
number of sources that may drive a system away from 
it. There are some mathematical approaches to study 
the effects of off-criticality on SLE based upon probabil- 
ity theory and conformal field theory 0, , which mainly 
pursue mathematical consistency in the continuum limit, 
but general feature of the off-critical driving function has 



not been known yet. In this Letter, we report our re- 
sults of numerical simulations to study how the Loewner 
driving function deviates from the ideal Brownian mo- 
tion when the system departs from the criticality in the 
case of percolation clusters. 

Loewner evolution. Let us start by reviewing basic 
elements of Loewner evolution briefly. Consider a non- 
intersecting continuous curve 7 which starts from the ori- 
gin and extends toward infinity in the upper half plane 
H. We parametrize 7 by t > 0, and denote a point on 
7 as 7t with 70 = 0. A part of 7 between 74^ and 74^ is 
represented by ^[t^^t2\- It is known that, for given 7[o.f], 
there exists a unique conformal map 5(4 : H \ 7(0, t] ^ HI 
that satisfies the condition 



gt{z) ^ z + at/z + 0{\z\-'^) 



(1) 
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with at > (Fig. [T]). The parameter t, which we will call 
time, is now defined by t :— at/2. The driving function 
Ut is defined by the image of 74 by the map gt- Ut := 
lim^^-yj gt{z). Then, it can be shown that gt{z) satisfies 
the Loewner evolution, dtgt{z) = 2/{gt{z) — Ut)- 

Note that gt{z) represents the complex electrostatic 
potential for the equipotential boundary of 7[o,t] ^nd the 
real axis [9|, then one can see that the time at — 2t is 
equal to p/2tt, where p is the 2-d dipole moment induced 
by 7[o,t] on a flat electrode. In this electro-static picture, 
Ut is just given by the charge induced by 7[o.t] along the 
right side of 7[o,t] and the positive part of the real axis in 
the unit system with Eq — I (Fig. [1]) . 

Schramm has shown that Ut becomes a Brownian mo- 
tion if the curve 7 is a conformally invariant random 
curve with the domain Markov property [l| . This means 
that any properties of the curve such as the fractal di- 
mension are determined solely by the diffusion constant 
of the Browninan motion. 

Simulations. We performed numerical simulations to 
obtain Ut for cluster boundaries of the site percolation 
on the triangular lattice with the occupation probability 
p < Pc — 0.5. Percolation clusters are generated in the 
rectangular system. To ensure the cluster boundaries to 
extend from the origin in the upper half plane, the left 
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FIG. 1: (Color online) A schematic diagram of conformal 
mapping, gt maps H \ 7(o,t] to H with 7t to Ut on the real 
axis. The deformed grid in the left panel can be regarded as 
equipotential and electric force lines, and is mapped to the 
straight grid in the right panel. 





FIG. 2: (Color online) Percolation on the triangular lattice. 
The left panel shows a whole system with left (right) half of 
the system boundary being occupied (unoccupied) . The right 
panel shows only the cluster boundary that starts from the 
origin. 



FIG. 3: (Color online) Examples of cluster boundaries (a- 
c) and their driving functions (d) for p = 0.5, 0.49, and 0.48. 
The boundaries are shown up to A'^ = 20000 steps. The lattice 
spacing is set to be unity. The correlation lengths are ^ = 400 
and 150 for p = 0.49 and 0.48, respectively. 



(right) half of the boundary sites are set to be occupied 
(unoccupied), and the boundaries are defined along the 
edges of the dual lattice, i.e., the honeycomb lattice 
(Fig. [2]). We use only the part of cluster boundaries that 
never touches the peripheries except for the bottom side. 
For each boundary curve, we numerically generate the 
conformal map gt{z) and compute the driving function 
Ut using the zipper algorithm 

Fig. [3l^a-c) shows typical examples of boundaries 7 
up to 20000 steps at p = 0.5,0.49, and 0.48, and Fig. 
[3]Jd) shows corresponding driving functions Ut p^ . For 
p = 0.5, the boundary extends into EI indefinitely around 
the origin once it hits an infinite cluster. For p < pc, 
the boundary tends to extend toward left in larger scale, 
but it can hardly be distinguished from that at p = pc 
within the scale of correlation length ^ with a finite lattice 
constant. Accordingly, Ut wander around the origin at 
p = Pc whereas it drifts toward —00 for large t for p < pc- 
Note that the time does not increases uniformly with the 
step along a boundary because the time is proportional 
to the dipole moment induced by 7[o,t]- 

The averaged behavior of Ut is shown for several values 
of p in the inset of Fig. IH where one can see {Ut) appears 
to drift at a constant velocity Vd- In Fig. [H the drift 
velocity Vd averaged up to t = 21000 over 4000 samples 
is plotted against pc — p in the logarithmic scale; Vd is 
shown to be proportional to (pc — p) for {pc — p) ^ 0.01, 

Vd-^{Pc~p), (2) 



whereas it behaves as 

Vd ^ {Pc - pT (3) 

with V « 4/3 for Pc - P > 0.01. 

The time dependence of the variance Var[C/t] is shown 
in the logarithmic scale in Fig. [5l As one can see in 
the plot for p = 0.42 in Fig. El^a), there are two time 
regimes: the superdiffusion regime and the normal diffu- 
sion regime, 

with the exponent a > 1 and the crossover time tc- 
The plots for several values of p in Fig. ^h) shows the 
tendency that the exponent a decreases to 1 and the 
crossover time tc increases as p pc- The estimated val- 
ues for tc are plotted against {pc — p) in the logarithmic 
scale by filled circles in Fig.lGja). 

Discussions. Our results can be understood as in the 
following. First, we consider the curve of p = Pc- In 
this case, 74 explores along a critical boundary, which 
extends indefinitely without bias. Consider a part of the 
curve 7[o,t] ■ This may look like a blob, whose size we 
denote by It- The time t that corresponds to the blob 
can be estimated as 

t-ll (5) 
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FIG. 4: (Color online) Drift velocity Vd vs. (pc—p). The drift 
velocities of the driving function averaged up to t = 21000 
are plotted in the logarithmic scale as a function of {pc — p) ■ 
There is a crossover around (pc — p) ~ 0.01 between the two 
regimes. The solid (dashed) line with the slope 1 (4/3) shows 
the limiting behavior in the small (large) Pc — p regime. The 
inset shows t vs. (Ut) for p = 0.5,0.4975,0.495,0.4925 and 
0.49 (from right to left). Each line represents average behavior 
over 4000 samples. 
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FIG. 5: (Color online) Var[?7t] vs. t in the logarithmic scale, 
(a) The plot for p — 0.42. The behavior changes from su- 
perdiffusive to diffusive around the crossover time tc indicated 
by the arrow. The behavior for p — pc, Va,r[Ut]~6t, is plotted 
for comparison, (b) The plots for several values of p. The 
crossover time tc becomes larger as we approach p^. 



using the electro-static analogy; The time t is propor- 
tional to the dipole moment p induced by j[o,t] attached 
to a flat electrode, and both the induced charge and the 
charge displacement are of order of It ■ On the other hand, 
Ut may be approximated to be Xt := Re7t because Ut is 
determined from the number of electric force lines land- 



FIG. 6: (Color online) (a)tc, ^, and vs. (pc^p) in the loga- 
rithmic scale. The horizontal solid line at t — 21000 indicates 
the simulation range, (b) —{xt) vs. (pc — p) for t = At, 2Af, 
and 3At with At = 5200, in the logarithmic scale. Each data 
point represents average over 3000 samples. 



ing on the right side of 7[o^f] and the positive part of real 
axis. Since xt wanders along the blob of size It, one can 
see that Ut undergoes normal diffusion without drift: 

Ut^lt- Vt. (6) 



For p < Pc, the curve 7 is the boundary of finite off- 
critical clusters that are connected by the bottom. Since 
typical size of each cluster is the percolation correlation 
length ^, the curve 7 looks like a chain of the blobs of 
size ^ (see for example Fig. [Sjc)). Within each blob, 7 
is almost like a critical boundary but with a bias toward 
the left; the strength of the bias is proportional to Pc — p- 

There should be two time regimes: the short time 
regime where 74 is still in the first blob, and the long time 
regime where 74 is traveling over blobs. The crossover 
time tc between the two regimes is the time that 7t goes 
over the first blob of size ^, therefore, it is estimated as 
tc ^ as in the case of critical blob. Within the short 
time regime, 74 explores in the first blob that looks al- 
most like a critical blob in the smaller length scale than 
^ but with a small bias toward left. We assume the 
scaling behavior (xt) ~ ~£,{t/tc)^ with an exponent p. 
If we determine the exponent so that the effect of bias 
should be proportional to {pc — p) for a fixed t, we obtain 
f3 = (v+\)j1v = 7/8, using the critical exponent v = 4/3 
for the correlation length of the 2-d percolation. On the 
other hand, the behavior of (xt) in the long time regime 
may be obtained by considering the situation where 74 
is at n-th blob, i.e., (xt) ~ n^. The time corresponding 
to this is t ~ n^^ because the dipole moment induced at 
each blob is of order of with which we can express 
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{xt) as t/^. With all these argument, we finally obtain 
i^t) - { (7) 

[-{Pc-prt it>u) 

with tc ^ (pc-p)^^". 

Within the approximation Ut ~ Xt, this is consistent 
with our results for the drift velocity Vd in Fig. 21 The 
slope we obtained for Pc — P }t 0.01 is very close to v ~ 
4/3. The linear dependence for Pc ~ P ^ 0.01 should 
correspond to the short time behavior. The crossover 
around pc—p ~ 0.01 is due to the fact that our simulation 
time length is not long enough in comparison with tc- 
If we can simulate longer time, the crossover value of 
Pc — P will be smaller. The expected weak non-linear t- 
dependence in Xt is difficult to distinguish from the linear 
behavior numerically; Actually, one might notice slight 
convexity of the plots in the inset of Fig. [H 

In order to check further consistency of the data with 
our interpretation, we plot tc for Var[f7t], and 
against Pc — P in Fig- El^a). The horizontal solid line 
indicates the time range of our simulations. The plotted 
range of the crossover time (filled circles) in Var[J7t] is 
too narrow to determine its behavior in a reliable way, 
but it seems consistent with that of (open circles). 
The value of Pc — P where the extrapolated tc of Var[C/(] 
reaches the simulation time is about 0.04 (arrow), which 



is larger than the crossover point in Vd by the factor 4; 
This discrepancy may come from the difference between 
Vd and Var[J7t], or simply due to the uncertainty in the 
estimates. In Fig. [Hb), — (xt) for several values of t 
are plotted against (pc — p). From Eq. ([7]), we expect 
linear dependence on {pc — p) for [pc — p) ^ t^^/^ and 
(Pc — p)^/^-dependence beyond that; The data seem to 
be consistent. 

With these observations, we interpret our results as 
follows. For an off-critical percolation cluster boundary, 
the driving function Ut undergoes a random walk with 
drift. The drift velocity Vd appears to be proportional 
to (pc — p) when [pc — p) ^ 0.01, but it is because time 
lengths of our simulations are finite; With the approxi- 
mation Ut Xt, the averaged driving function should be 
given by Eq. ([7]). Within the crossover time, the random 
walk is superdiffusive with exponent a > 1, which de- 
creases toward 1 as p — > pc- Beyond the crossover time, 
the fluctuation around the drift motion is normal diffu- 
sion with larger diffusion constant than that at p = Pc. 

Finally, let us discuss the scaling limit where the lattice 
constant a — > with keeping the correlation length f 
constant. Our argument to derive Eq.(7) holds for any f 
and a as long as ^ 3> a, thus if we scale the variables as 
it = Xt/£,, Ut = C^t/^, and i = t/^^, then we obtain the 
corresponding equation for the scaled variables, i.e. the 
same as Eq.(7) but without the factor of (pc — p). 
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